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Abstract 

The aim of this work is to analyze the entropy, entropy flux and entropy rate of granular materials 
within the frameworks of the Boltzmann equation and continuum thermodynamics. It is shown that 
the entropy inequality for a granular gas that follows from the Boltzmann equation differs from the 
one of a simple fluid due to the presence of a term which can be identified as the entropy density 
rate. From the knowledge of a non-equilibrium distribution function - valid for for processes closed to 
equilibrium - it is obtained that the entropy density rate is proportional to the internal energy density 
rate divided by the temperature, while the entropy flux is equal to the heat flux vector divided by the 
temperature. A thermodynamic theory of a granular material is also developed whose objective is 
the determination of the basic fields of mass density, momentum density and internal energy density. 
The constitutive laws are restricted by the principle of material frame indifference and by the entropy 
principle. Through the exploitation of the entropy principle with Lagrange multipliers, it is shown 
that the results obtained from the kinetic theory for granular gases concerning the entropy density 
rate and entropy fiux are valid in general for processes close to equilibrium of granular materials, 
where linearized constitutive equations hold. 

PACS: 45.70.-n; 51.10.+y; 05.70.Ln 
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1 Introduction 

In the last two decades a great number of papers has appeared in the literature concerning the theory and 
applications of granular gases. The determination of the equilibrium and non-equilibrium distribution 
functions from the Boltzmann equation was investigate among others by Lun et al. [T], Jenkins and 
Richman [5J |3] , Goldshtein and Shapiro [i] , Brey et al. [5] and a more complete list of the works can 
be found in the books [3 El |H] ■ Whereas for a rarefied gas of elastic particles the collisions conserve the 
mechanical energy and the gas relaxes to an equilibrium state described by a Maxwellian distribution 
function, the inelastic collisions between the particles of a granular gas transform the translational kinetic 
energy into heat and there exists no equilibrium state which is characterized by a Maxwellian distribution 
function. Several remarkable properties of granular gases are reported in the literature, among others, 
the decay of the temperature of a granular gas (the so-called Haff's law), nonhomogeneous structure 
formation (cluster formation), shear instabilities (shock wave formation), anomalous diffusion, etc. 

While the transport properties of granular gases are almost well understood, there is a lack in the 
literature, to the best of our knowledge, concerning the discussion of the entropy, entropy flux and entropy 
rate of a granular material and the aim of this work is to analyze these subjects within the framework of 
the Boltzmann equation and from the viewpoint of continuum thermodynamics. 

In Section 2 we obtain the balance equation for the entropy density of a granular gas from the 
Boltzmann equation, by showing that there exists a positive semi-definite quantity which can be identified 
with the production rate of the entropy density. The so-called entropy inequality for a granular gas differs 
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from the one of a simple fluid due to the presence of a term which can be interpreted as the entropy 
density rate. Due to the energy loss of the gas particles an entropy density rate must be present in order 
to preserve the positiveness of the production rate of the entropy density. From the knowledge of the 
non-equilibrium distribution function - which is valid for processes close to equilibrium - the constitutive 
equations for the pressure tensor, heat flux vector, internal energy density rate (which is related to the 
cooling rate of a granular gas), entropy flux and entropy density rate are calculated. It is shown that 
the entropy density rate is proportional to the internal energy density rate divided by the temperature, 
while the entropy flux is equal to the heat flux vector divided by the temperature. 

The results obtained in Section 2 are used in Section 3 as a motivation for the development of a 
thermodynamic theory for granular materials. The objective of such theory is the determination of 
the basic fields of mass density, momentum density and internal energy density. The necessary field 
equations are based upon the balance equations of the basic fields which are closed by constitutive laws. 
The principle of material frame indifference and the entropy principle are used in order to restrict the 
generality of the constitutive equations for the internal energy density, pressure tensor, heat flux vector, 
internal energy density rate, entropy density, entropy flux and entropy density rate. From the exploitation 
of the entropy principle with Lagrange multipliers we show that for processes close to equilibrium, the 
relationships between the entropy density rate and the internal energy density rate and between the 
entropy flux and the heat flux vector found in Section 2 by using a restricted non-equilibrium distribution 
function are valid in general for granular materials within a linearized theory. 



2 Kinetic theory of a rarefied granular gas 

In this section we review briefly the main features of a granular gas which is described by the Boltzmann 
equation and obtain the balance equation for its entropy density. We consider a rarcfled granular gas 
free of external body forces where only binary collisions between the particles are taken into account. 
Let us denote by m and d the mass and the diameter of a spherical particle and by (c, Ci) and (c',Cj) 
the velocities of two particles before and after collision, respectively. If g = Ci — c and g' = c'^ — c! 
represent the relative velocities before and after collision, the inelastic collisions will be characterized by 
the relationship (g' • k) = — a(g • k) where < a < 1 denotes a normal restitution coefficient and k is the 
unit collision vector which joins the centers of the two colliding spheres pointing from the center of the 
particle denoted by 1 to the center of the other particle without index. 
The velocities before and after collision are related by the equations 

c' = c + ^— (g-k)k, c'l-ci ^(g-k)k. (1) 

From the above equations it follows the relationships which connect the relative velocities and their 
modulus, namely, 

g' = g-(l + a)(g-k)k, g'^ = g' -{\-c?){^-Vf. (2) 
Furthermore, the variation of the kinetic energy in terms of the pre- and post-collisional velocities, reads 

_c'2 + _cf _ _ = _(a2 - l)(g . k)2. (3) 

When a = 1, which is the case of elastic collisions, it follows the conservation of the kinetic energy. 

A direct encounter is characterized by the pre-coUisional velocities (c,Ci), by the post-collisional 
velocities (c'jC'^) and by the collision vector k. For a restitution encounter the pre- and post-collisional 
velocities are denoted by (c*, c^^) and (c, Ci). respectively, and the collision vector by k* = — k. It is easy 
to verify that the relationship g • k* = — Q!(g* • k*) = — (g • k) holds. 

In a restitution encounter the pre-collisional velocities (c*,c*) are related to the post-collisional ve- 
locities (c,Ci) by the equations 

c = c* + ^— (g' •k*)k*, ci = c^ ^^(g*.k )k*. (4) 
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By using the relationships k* = — k and (g • k) = — a(g* • k) that characterize a restitution colhsion, the 
above equations can be written as 

c*=c+i^(g.k)k, =ci-i^(g.k)k. (5) 

For the determination of the Bohzmann equation we have to know the transformation of the volume 
elements dcj rfc* = \J\dcidc where |J| is the modulus of the Jacobian of the transformation. By a 
straightforward calculation \J\ = 1/a and it follows that 

(g* • k*) rfc* dct = ^(g • k) dcdci. (6) 

From the expression ^ one may infer that the Boltzmann equation for granular gases without external 
forces is given by 

where x is the radial equilibrium distribution, which takes into account the correlation between the 
positions of the colliding particles. 

The multiplication of the Boltzmann equation ([7|) by an arbitrary function ip(x,c,t) leads to the 
so-called transfer equation 
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dtp dtp 
dt ' dxi 



fdc 



J ^(x,c,t) (^^f*f*-f,f^ d2(g.k)dkrfcidc 
= I y"[7^(x,c;,0+V'(x,c',i)-7^(x,Ci,i)-^(x,c,0]/i./d'(g-k)dkdcidc. (8) 

The second equality above was obtained by taking into account the relationship ([5]), by renaming the 
the pre-coUisional velocities (c*,c*) as (c,Ci) and the post-collisional velocities (c, Ci) as (c',c'j^) and by 
using the symmetry properties of the collision term when the two molecules are interchanged. 

A macroscopic state of the granular gas is characterized by the fields of mass density g, hydrodynamic 
velocity and internal energy density ge defined by 



e = J mf dc, gvi = J maf dc, ^ J 

where Ci — Vi is the so-called peculiar velocity. 

The balance equations for the fields © are obtained by choosing tp equal to m, mci and toC^/2 in 
the transfer equation ([8]), yielding 

■ ^^"'-0, (10) 



dt dxi 
dgvj _^ d{gviVj +Pij) 
dt dxj 



= 0, (11) 



dge d {gev^ + g,) dvi . ^ 

^ + d^. = (^2) 

Above, Pij and qi denote the pressure tensor and the heat flux vector, respectively, which are defined by 

P^J = J mac, fdc, q, = J "^C^CJdc. (13) 

Furthermore, C is the internal energy density rate of the granular gas. From the right-hand side of ([8]) 
and the relationship ([3]) it follows that 

C = / /i/(g • k)3dkdci dc. (14) 
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The so-called cooling rate (* is related to the internal energy density rate C C* ^ 2mC,/3kTg. Note 
that the internal energy density rate vanishes for elastic collisions, i.e., for a = 1. 

By applying the Chapman-Enskog method to determine the non-equilibrium distribution function 
from the Boltzmann equation (O it follows [S] [Hj 

^ g/rnxi^_^[ 16(l-a)(1^2a^) f 15 5PC^ B^C^\ 

■' m\2TTkTJ 1 81-17a + 30a2(l-Q) V 8 2 2 J 



Above, k is Boltzmann's constant and the temperature is related to the specific internal energy through 
T = 2me/3fc. Moreover, the coefficients 71, 72 and 73 read 



15m / m —15m 1 / m \3 

Xc|2gT(9 + 7a)(a + l) V " 2xd^g{li -a)(a + l)\n \kf) ' 



300m(l - a) jm 
X^^Q^i'd + 7a){a + 1)(19 - 3a) \l M^' 



We may observe that the non-equilibrium distribution function (fT7|) does not reduces to a Maxwellian 
distribution function in the absence of the spatial gradients due to the presence of the underlined term. 
However, this term vanishes in the case of elastic collisions, i.e., for a = 1 . 

The insertion of the distribution function ([T5|) into the definitions of the pressure tensor pij and heat 
flux vector qi - given by (|13p - and subsequent integration of the resulting equations lead to the linearized 
constitutive equations 

= 9-TS,j - 2m ( ^ - -g-S,, 1 , q.^ = -no—-^i — . (18) 

Above, the parentheses around the indexes denote the symmetric part of a tensor. Furthermore, the 
coefficients of shear viscosity thermal conductivity and ki are given by 



15 mkT 75 k mkT 

^ ~ 2xd2(13-a)(a + l)V tt ' ~ 2xd^(9 + 7a)(a + 1) m V tt ' 

750(1 -a) kT fkT , , 

Ki = — ^ ^ ^ \ . 20 

Xd^(9 + 7a)(a + l)(19-3a) g \l mir 



The linearized expression for the internal energy density rate Ci obtained from the insertion of ([15 
into (jl4p and subsequent integration of the resulting equation, yields 



Co = 2xfcT^d^/^(l-a2) 
m^ V m 



1 



3(l-a)(l-2a2) 
81 - 17a + 30^2(1 - a) 



(21) 



The above expression for the internal energy density rate is valid only in a linearized theory where 
products of gradients are neglected. 

Now we may choose ip{x,c,t) — — fcln/ into the transfer equation ([5]) and obtain 



d f d f Yd k f f f 

- J i-k\nf)fdc+— j {v, + a)i~k\nf)fdc+^ J ln^/i/'(g.k)dkdcidc 

= ^ / fj ( A7 " ^'^^^ ■ (22) 



Equation ([^^ can be written as 

dgs d{gsv^ + ip,) 
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which can be identified with tlie balance equation for tlie entropy density gs, where ipi is its entropy flux, 
c; denotes the entropy density rate and a its production rate. The expressions for these quantities are 
given by 



gs = -fc//ln/dc, c^^^ In ij^fi fig. k)dkdc,dc, (24) 



ip, = -kjcjlnfdc, a=^llnj^(^j^-?jhf{g-k)dkdc,dc. (25) 

On the basis of the relationship (x — l)lna; > 0, which is valid for all x > 0, we may infer that the 
production rate of the entropy density is a non- negative quantity, i.e., a > 0. 

From the insertion of the distribution function (|T5]) into the definition entropy flux (|25)) i and the 
integration of the resulting equation leads to the following linearized expression 

V'. = |. (26) 

Hence, in a linearized theory the entropy flux of a granular gas is equal to the heat flux vector divided 
by the temperature. 

The determination of the entropy density rate (j24p 9 and its production rate (125^ 9 is more involved. 
First, we may write the entropy density rate as 

, = ^ 1 In jjlifis ■ k)dkdc,dc = ^ / In jjfins' ■ k')dkdc[dc' 

= ^ / In j^hfis ■ k)dkdcidc, (27) 

by using the relationship (g-k)(icrfci = -^{g' -k') dc' dc[ and by renaming the velocities (c',c'j^) i — > (c, Ci) 
(c, Ci) I — > (c*, cl). In a linearized theory 

'» # - 2if '^f + - - ^ iif <' - "''fe- ■ ^'>^ = iif (^) - 

The insertion of (|28p together with the distribution function (|15p and integration of the resulting equation 
leads to 

showing that the entropy density rate in a linearized theory is proportional to the internal energy density 
rate. 

By following the same methodology, it is easy to obtain that 

In ^ = -^(C? + C'^ Cl -C^) = ^il- a2)(g . k)2^ (30) 
so that production rate of the entropy density in a linearized theory reduces to 



« = -r^, (29) 



Note that for restitution coefficients a < 1, the internal energy density rate ([2T|). the entropy density 
rate ([29]) and the production rate of the entropy density ([3T|) are always positive semi-definite quantities, 
i.e., Co ^ 0, ijo > and uo > 0, with the equal sign valid in the elastic limit, i.e., for a = 1. For restitution 
coefficients close to unit (a i — > 1) the internal energy density rate (PT|) becomes 



Co = AxkT^A\^{l -a) + 0{{1 - af), (32) 
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and the entropy density rate ([29|) and its production rate ([31]) reduce to 



^o = ^ + 0((l ao = 0((l-a)2). (33) 

Hence for quasi-elastic restitution coefficients the entropy density rate is equal to the internal energy 
density rate divided by the temperature, while the production rate of the entropy density becomes a 
term of second order in (1 — a). 



3 Thermodynamics of granular materials 



The main objective of a thermodynamic theory of granular materials is the determination of the fields 
of mass density p(x, i), velocity Wi(x, t) and temperature T(x, i) in all points of the material x at time 
t. The knowledge of these fields are based on the balance equations of mass density pop . momentum 
density (fTT|) and internal energy density ([T^ . 

However, the system of equations pUj) through ([T^ is not closed for the determination of the basic 
fields g, Vi and T, since it is necessary to express the constitutive quantities pij, qi, e and in terms of the 
basic fields. The generic form of the constitutive equations which obey the principle of material frame 
indifference and refer to a viscous and heat-conducting granular material is given by 



{s,qi,PijX} ^ J^i Q,T, 



dT 



dQ_ 

dx,. 



dx 



(34) 



The solutions of the balance equations ([T(I|) through ([T^ for the basic fields which takes into account the 
constitutive equations (j34p arc called thermodynamic processes. 

The constitutive equations are restricted also by the entropy principle which states that the entropy 
inequality ()23p must hold for every thermodynamic process. Furthermore, the specific entropy density s, 
the entropy flux Lpi and the entropy density rate ^ are considered constitutive quantities whose dependence 
is the same as the ones given in ([M)) . namely. 



dQ_ dv^ 

OXi OXi OXj) 



(35) 



In this work we are interested in a linearized theory, so that we may write the following linearized 
constitutive equations 

dxi ' 
dxi ' 



c 



£o(f,T)+£i(f>,T) 

Co(^?,^) + Cl(^^,T) 



dxi 
dvj 
dxi ' 



- so{g,T) + si{g,T)- 



-Ko{g,T) 



dT 



Pro 



dx. 

p{q,T) - i]{g,T) 



ni{g,T) 

dVr 
dXr 



dg 



dxi' 
% - 2fi{g,T) 



= -<poig,T) 



dT 



dvi. 



dx 



J) 



dxi 

1 dVr r 
3d^r '' 



dg_ 

dxi 



(36) 
(37) 
(38) 
(39) 



Note that we are dealing with a granular material where the specific internal energy density e, the 
specific entropy density s, the internal energy density rate C and the entropy density rate ^ may depend 
on the mass density, temperature and on the divergence of the velocity. Furthermore, in the constitutive 
equation for the pressure tensor, rj is the coefficient of bulk viscosity. This term does not appear in (|18|)i 
due to the fact that in Section 2 the granular material refers to a rarefied monatomic gas. 

The exploitation of the entropy inequality proceeds by using the method of Lagrange multipliers [SI [IH] 
and by imposing that the inequality 



dgs 
'dt 



d {gsvi + ipi) 



dxi 
dge 
~dt 



d{gevi + qi) 
dxi 



dgvi 
dx, 
dv. 



-A, 



/ dgui d{gViVj + Pij) 



V dt 

>0, 



dxi 



(40) 
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must hold for all thermodynamic processes. Furthermore, the Lagrange multipliers A, Ai and A are 
considered functions of 

{A,A,,\} = T[g,T,—,—,- 



dxi dxi ' dx 



3) 



(41) 



The insertion of the constitutive equations p6p through (|39p into the entropy inequality (^U)) leads to 
an inequality which is linear in the derivatives 



dvi 
'dt 



dg 
dt' 



dT 



dvk 



k.k 



dt 



dvk,k 
dxi 



dxi 



d0,j 
dxi 



(42) 



where the comma indicates the differentiation with respect to the spatial coordinates. The resulting 
inequality must hold for all values of the quantities (|42|) . so that the coefficients of these derivatives must 
vanish and we get 



ds 



/9s ^de 
^ \ dg dg 



A = 0, 
= 0, 



Furthermore, there remains the following residual inequality 



dg 



ds de 
df^df"' 



dg, 



dg dg 



dg 

dx, 



dtpi _ dqi 
dT dT 



dT 

dxi 



' dxi 



dv,. 



— + <; - AC - Ag^ - Xp^ — > 0. 



(43) 
(44) 

(45) 



Now wc proceed to analyze the conditions (^51) and pij) . First the insertion of the constitutive 
equations (|38]) for the heat flux vector and entropy flux into the equations (|44)) 9^ leads to 

(46) 



ipi = Ag,;, where A = X{g,T). 
Next from equations (03]) 3, (j44l)i and (|46])2 it follows that 



djs - As 
dvk.k 



0, 



djs - Xs) 
df 



dX 
'df' 



(47) 



From the above equations we may conclude that the specific internal energy density e and the specific 
entropy density s do not depend on the divergence of the velocity and in a linearized theory the constitutive 
equations ([551 reduce to 

e^eo{g,T), s = so(e,T). (48) 



These results are the same as those found for a simple fluid [TOj [TT] . 

The next step is to build the differential of the specific entropy density s, which by the use of the 
relationships (|i5))9^. yields 



ds 

ds = —dT - 
dT 



ds_ 

dg 



dg = X \ de 



A_ 

'qX 



dg 



The comparison of the above equation with the Gibbs equation of thermodynamics 



ds 



1 

f 



de 



jdg 



leads to the identification of the Lagrange multipliers 



A 



1 



A 



'Tg 



(49) 



(50) 



(51) 



From (|46Pi we get that in a linearized theory the entropy flux is equal to the heat flux vector divided 
by the temperature, i.e., (pi — qi/T, which is the same result as the one of the kinetic theory of rarefied 
granular gases. 



7 



For the exploitation of the residual inequality we insert the constitutive equations ([36]) through p9|) 
and the relationships ([51]) into the residual inequality (|45|) and obtain 



if dT dg\ dT dvi 1 / dvi\ rj dvi dvj ^fJ. dv^i dv^i 

T^[''°d^,+'''d^Jd^,+'° + ''d^rT[^' + '''d^J + T ^^^^ 

Above, we have introduced the velocity gradient deviator defined by 

dxj) dxj) 3 dxr 

As it was pointed out in the last section, the production rate of the entropy density (PT|) does not 
vanish unless a = 1, which correspond to the elastic limit. In the elastic limit the internal energy density 
rate and the entropy density rate vanish and the granular material reduces to a single fluid. Hence, in 
a state where the spatial gradients vanish we get from the residual inequality ((52)) that <ro > Co/T, with 
the equality sign valid only when the production rate of the entropy density vanishes, i.e., in the elastic 
limit. As was expected, this conclusion is the same as the one found from a kinetic theory of rarefied 
granular gases (see ([29| . ([3T|) and their approximations ([33]) when a i — > 1). Furthermore, the residual 
inequality (j52p is linear in the divergence of the velocity and the coefficient of this term must vanish in 
order to preserve the inequality and it follows that <;i = Ci/T- Hence, as in the kinetic theory of rarefied 
granular gases, if the production rate of the entropy density could be consider as a term of small order, 
the entropy density rate in a linearized theory is equal to the internal energy density rate divided by the 
temperature, i.e., = C,/T when cr w 0. However, note that in the kinetic theory of gases this conclusion 
is valid only when the restitution coefficient is close to one, which corresponds to a quasi-elastic collision. 
It is worth to call attention to the fact that the entropy density rate and the internal energy density rate 
for granular materials in a linearized theory depend also on the divergence of the velocity. 

It is not possible to extract more information on the coefficients of shear and bulk viscosities and 
thermal conductivity from the residual inequality ()52p . since we have two scalar constitutive quantities 
- namely the internal energy density rate C and the entropy density rate <; ~ and we have considered 
linear representations for both. By taking into account quadratic terms for these two quantities, the 
exploitation of the residual inequality leads to results that are not remarkable to discuss here. 

It is interesting to analyze the case of a simple fluid where the entropy density rate and the internal 
energy density rate are absent. In this case we conclude from the residual inequality that the coefficients 

kq ki/2' 
Ki/2 

positive semi-definite matrix. This last result implies that the coefficient K2 must vanish for a simple 
fluid and Fourier's law reduce to q = —kq\7T. This conclusion is the same as that obtained by Liu [TT] 
by analyzing a simple fluid. 



of shear and bulk viscosities are non- negative, i.e., /i > and rj > 0^ and that 



4 Final Remarks 

It is important to call attention to the fact that the relationships obtained from the thermodynamic theory 
in Section 3 are not restricted to rarefied monatomic granular gases as in Section 2 but also valid for 
granular materials in general, i.e., monatomic, polyatomic and real gases. Indeed, in the thermodynamic 
theory there is no dependence of the constitutive equations on the restitution coefficient, the pressure 
tensor has a non- vanishing bulk viscosity term and the entropy density rate as well as the internal energy 
density rate depend on the divergence of the velocity. The dependence of pressure tensor, the entropy 
density rate and the internal energy density rate on the divergence of the velocity could be achieved within 
the framework of the kinetic theory of polyatomic or dense granular gases, since their non-equilibrium 
distribution functions depend on the divergence of the velocity. 
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